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Abstract 

We investigated the physical mechanism for the pattern transition from 
square lattice to stripes, which appears in vertically oscillating granular lay- 
ers. We present a continuum model to show that the transition depends on 
the competition between inertial force and local saturation of transport. By 
introducing multiple free-flight times, this model further enables us to analyze 

the formation of superlattices as well as hexagonal lattices 
PACS numbers: 47.54.+r, 45.70.QJ, 45.70.-n 



Typeset using REVTgX 



* e-mail : childend@complex. kaist . ac . kr 
TEL:+49-351-871-1225 FAX: +49-351-871-1999 

1 



Vertically oscillated granular layers exhibit various interesting patterns such as sub- 
harmonic stripes, squares, hexagons as well as localized structures called oscillons 
The control parameters are the driving frequency / and the dimensionless acceleration 
T = 4n 2 f 2 A/g, where g is the acceleration of gravity, and A is the amplitude of the vibrating 
surface. Increasing T with fixed /, flat states, subharmonic squares or stripes depending on 
/, subharmonic hexagons and flat states appear in sequence. 

These granular patterns have been studied using phenomenological models []8| |3~5T and 



molecular dynamics simulations [ 16fl . The studies show that transitions between granular 



patterns as F increases are related to the period doubling bifurcation P,|12] . But the bifurca- 
tion as a function of /, that is the transition between squares and stripes, has not been well 



understood yet. There have been many models showing both squares and stripes |8| JTQ]T^JT6 
and it was suggested that the effective viscosity and the horizontal mobility resolve the tran- 
sition [plHH. Experiments also show that squares arise at large lateral transport while stripes 
form at small lateral transport |5|]. However, there has been no explicit model for the square 
to stripe transition upon increasing frequency. 

In this paper, we propose a continuum model which describes the transition between 
squares and stripes. We incorporate the lateral transfer and the local saturation mechanism 
into our model. The study shows that the ratio of these two competing processes depends 
on the frequency, which explains the dependence of the square-stripe transition on the drive 
frequency. In addition, considering multiple free-flight times, we study the hexagonal and 
superlattice patterns which are observed in the granular system with multiple frequency 
forcing M. 



The granular patterns are similar to Faraday waves in vibrated fluids , however some 
differences exist. 

(i) On each cycle, there is a free-flight time during which granular layer loses contact with 
a bottom plate. 

(ii) In general, macroscopic physical quantities for the granular material are not well de- 



fined. But Ref. shows that a continuum description of the granular layer is valid 
for the square-stripe transition. Molecular simulation results [[RJ also suggest that 
granular continuum equations can be used as far as pattern formation of vertically 
oscillated granular media is concerned. 

(iii) Granular material has an angle of repose, and shows both solid-like and liquid-like 
motions. 

Considering (i),(ii), a continuum equation is introduced with two steps: contact with 
bottom plate, and free- flight. It provides the two competing processes, diffusion during the 
contact and focusing effect during the free flight, which were suggested as mechanisms of 
subharmonic wave generation ||. A hysteresis effect due to (iii), which accounts for the 
formation of oscillon P, p!0| , p!2| -[T5|] , is not considered here since it appears to be unrelated to 



the square-stripe transition [[J. 

At t = 0, the layer collides with the bottom plate and the slope of granular layers begins 
to relax. From the collision, granular particles get kinetic energy and move randomly. The 
process can be described as diffusion, and the following equation is obtained for < t < t c 
during which granular layers have contact with the bottom plate 

Pt = DV 2 p. (1) 

We assume that the thickness h is a monotonic increasing function of the local density p, 
where p is a mass of granular material per unit area, and choose p as an order parameter. 

At t = t c , the layer loses contact with the bottom plate. During the free flight, the lateral 
transfer of grains induce the focusing of the sands, so that subharmonic excitation of the 
granular wave occurs. From continuity, the flux j and the velocities v of grains right after 
takeoff from the bottom plate are given by 

j(x,t c ) = -DX7p(x,t c ) 

(2) 

v(x,t c ) = j(x,t c )/p(x,t c ) = -DVp(x,t c )/p(x,t c ). 



This is indeed what was done by Cerda et. al. ||, in which saturation of the flux is not 
considered, and only square patterns arise. Hence one may expect that the saturation is 
related to the formation of stripe patterns. Assuming that granular layers behave like fluid 
during the free flight, we adopt Navier-Stokes equations in 2 dimensions as a model of the 
the local saturation of flux as well as the transfer of particles and momentum. 

During the interval t c < t < T, where T = 1/f is the period of the driving force, we get 

p t + V • (pv) = 

(3) 

v t + v ■ Vv = —^Vp. 

Momentum is transferred through the transport of the granular particle, and inelastic 
collisions between the particles suppress the lateral motions. In general, these effects can be 
incorporated into the local stress tensor 7y. Assuming that the suppression of the lateral 
motion is most important, the local stress tensor is approximated by T^- = —p5ij. When we 
consider the collision between granular particles, the relevant variable is a relative velocity 
rather than a velocity itself. A dimensional consideration leads to the following form of the 
pressure 

p = — ap\ V • v\V ■ v. (4) 

The pressure is high at the accumulation points where V ■ v < , and low at the dispersion 
points where V • v > 0. The pressure difference produces the force to suppress the focusing 
of sands. The parameters of Eq. (0),(@), (§) and (|J) are D, f,a and t — tf/T — 1 — t c /T 
where tf is a free flight time. In a real system, r is determined mainly by T. D and a 
depend on layer depth, grain diameter, and restitution coefficient as well as T and /. We 
choose a = 0.015, D = 0.17, and system size I = 10. 

Linear stability analysis for the rest state shows that the amplitude of the mode with wave 
number k is amplified by the Floquet multiplier F(k) = (1 — Dk 2 rT) exp(— Dk 2 (l — r)T) on 
each cycle. At small r, we get only flat states. As r is increased beyond r c = 0.782, which 
is independent of D and T, minimum value of F(k) becomes less than -1 and subharmonic 
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standing waves arise. In this case, we obtain the dispersion relation, / = Dt(1 — r)k 2 . 
Note that the value of r c does not depend on D, f,a and I. The bifurcation diagram of Eq. 
([[]) and (|]) as a function of / and r is presented in Fig.|TJ; Fig.^| describes typical patterns. 
With increasing f with fixed a and r , one can see the square to stripe transition as in Fig.|2| 
(a),(b). 

Increasing a with fixed f and r also yields the square to stripe transition. In other words, 
the local saturation prefers stripe patterns. It is interesting to compare the results with the 
Faraday waves. In Faraday waves, squares are observed for small viscosity while stripes for 
large viscosity ||17|| . In this case viscosity represents the strength of local saturation which 
competes with the inertial force. 

To understand the relation between the square-stripe transition and the drive frequency, 
let us write the equations in dimensionless forms. Naturally, we take T = 1/f as the time 
unit. For the length unit, we choose J 1 D / f which reflects the wavelength. Inserting t = tj f, 
V = V ' J^JdJJ, v = ^Dfv, into Eq. (|1|), Eq. ([|) with Eq. ([|), and dropping the carets, we 
get for < t < 1 - r, 

Pt = V 2 p (5) 

and for 1 — r < t < 1, 

Pt + V ■ (pv) = 

(6) 

v t + v ■ W = f V(p|V • u|V • v) 

where R is defined as ^j. 

Now we have two dimensionless parameters. One is r — tf/T, the fraction of flight time 
for one period, and the other is R = ^ in Eq. (||), which corresponds to the ratio of the 
force causing the local saturation to the inertial force. Note that R is proportional to the 
frequency /. It means that increase in / induces the strong local saturation compared to the 
lateral movement. Qualitatively, this can be explained as follows. Lateral movements are 
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suppressed by the inelastic collisions which are reflected by pressure p in our model. When 
wavelengths are short, the velocity differences between the neighbor grains are large, then 
so is pressure p from Eq. (|J). Hence, for high frequencies, the collision is dominant and 
mobility is suppressed. 

Experimentally, the dispersion relation has a kink structure, i.e., it takes a different 
form below and above / = f c H5|JqJl6| . For almost all conditions, the square-stripe transition 
occurs for / < f c || . In other words, only the dispersion relation for / < f c is relevant as far 
as the square-stripe transition is concerned. In this regime, earlier experiments |]|| reported 
an exponent 2, however a recent experiment || yielded a different exponent value of 1.3. In 
contrast, in our model, the dispersion relation A ~ [D j f) 1 ^ 2 is obtained. This contradiction 
can be resolved by considering the frequency dependency of the diffusion coefficient D p|,|T5 



Considering the relation between viscous length and energy injection rates, for fixed T, D 
scales like l// 3 in the small frequency regime PJI8|. With this, A ~ (l// 2 ) i s obtained and 
the parameter R which measures the ratio of the saturation of mobility to the inertial force 
is proportional to f 4 . Also note that influences of T on D and a have not been considered. 
Dependence of the square to stripe transition frequency on r in Fig.|T], in contrast to the 
experimental result where the transition frequency depends only on the layer depth H, 
may be attributed to this. Nevertheless, for fixed T, the above conclusion that squares 
transit to stripes with increasing the drive frequency does not change. 

Next we consider the instability of stripe patterns. With an abrupt change of the driving 
frequency /, we observe a crossroll and a zigzag instability. Fig. |3](a) shows an example 
of the crossroll instability. The stripes are unstable to crossrolls on both the high and low 
k sides in our model. Fig. |3](b) illustrates the increase in wave vector k due to the zigzag 
instability. In experiments and molecular dynamics simulations ||19|| , the stripe patterns have 
the crossroll instability for decreasing k, and the skew-varicose instability for increasing k. 
The zigzag instability has not been observed in experiments. Refining the stress tensor and 
comparing the results with the instabilities in experiments will yield more insights about 
the interaction between granular particles. 
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As external forcing T is increased, the motions are repeated after two different flight times 
and contact times. This period doubling causes hexagonal patterns ||, and the mechanism 
is investigated theoretically fl2| . Recently, it has been reported that multiple frequency 
forcing on the granular layer yields superlattice patterns [[7]. Similar patterns have been 
studied in Faraday waves [f2t|, driven ferrofluid surfaces []21j] , nonlinear optics |2"2" |, Turing 



patterns |23| and Rayleigh-Benard convection [p4| . As in the single frequency forcing ||, 



the important variables determining patterns are expected to be the free-flight times and 
contact times. 

To study the hexagons and superlattice patterns, we introduce the multiple free-flight 
times tfi, tf 2 , t/ n and contact times t&, t c2 , ■■■■ l t cn . Let's apply Eq.([l]) during t cl and Eq.([3D 
during tf\. After that, we apply again Eq.(|l|) during t c2 , Eq.(|3]) during tf 2 and so on. The 
case with n = 2 reduces to the period doubling and hexagonal patterns are obtained. 

Linear stability analysis shows that the amplitude of the mode with wave number k 
is amplified by the factor F(k) = nr=i(l ~ Dk 2 tfi) exp(—Dk 2 t ci ) on each cycle. For the 
hexagonal patterns, as expected, F(k) gives a harmonic mode as in Fig. ^(a). However, in 
general, F(k) can yield instabilities with various wavelengths and, as a result, superlattice 
patterns can be formed. In Fig. |](b), the modes with k\ and k 2 are subharmonic and the 
mode q is harmonic. They satisfy the resonance condition k\ + k 2 + q = 0. Evidently, there 
is a large freedom of choice for the values of t C i and tfi. Analysing the characteristics of the 
superlattice patterns in these vast parameter spaces is a problem for future studies. 

In conclusion, based on the mass and momentum conservation law, we presented the 
continuum model to explain the transition between squares and stripes in vertically oscillated 
granular layers. Patterns are selected by two competing nonlinear interaction: the inertial 
force and the local saturation due to the inelastic collisions between granular particles. With 
increasing the drive frequency, the local saturation gets stronger compared to the inertial 
force effect, so that the square to stripe transition occurs. Introducing multiple free-flight 
times and contact times yields hexagonal patterns and superlattice patterns. 
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FIGURES 



FIG. 1. Bifurcation diagram as a function of / and t. a = 0.015, D = 0.17, I = 10 

FIG. 2. Square and stripe patterns. White corresponds to the large p. D = 0.17,/ = 10, 
a = 0.015, r = 0.85. (a) T = 1.2, (b) T = 0.7. 

FIG. 3. Instabilities of stripes, a = 0.015, D = 0.17,1 = 10. For t < 0, T = 0.55, and at t=0, 
T = 1/ f is changed abruptly, (a) T=0.85, Crossroll instability, (b) T=0.30, Zigzag instability. 

FIG. 4. Hexagon and superlattice patterns. Left : Hexagonal and superlattice pat- 
terns; Middle : Power spectra of the left images; Right : Floquet multiplier F(k) vs k. 
D = 0.17,1 = 10, a = 0.015 (a) t cl = 0.28, t c2 = 0.08, t fl = 0.84 and t f2 = 1.02. (b) t cl = 0.08, 
t c2 = 0.18,^3 = 0.31,t c4 = 0.03, t fl = 0.72,t /2 = 1.21,t /3 = 2.00 and t f4 = 0.15. 
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